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Abstract:  Molecular-level  modelling  and  simulations  of  the  high-pressure  volumetric  response 
and  irreversible  densification  of  a  prototypical  soda-lime  glass  are  first  employed.  The  molecular- 
simulation  results  obtained  were  next  used  to  modify  the  pressure  versus  degree-of-compression 
(the  negative  of  volumetric  strain)  and  yield  strength  versus  pressure  relations  in  order  to  account 
for  the  effects  of  irreversible  densification.  These  relations  are  next  used  to  upgrade  the  equation 
of  state  and  the  strength  constitutive  laws  of  an  existing  material  model  for  glass.  This  was  fol¬ 
lowed  by  a  set  of  transient  non-linear  dynamics  calculations  of  the  transverse  impact  of  a  glass 
test  plate  with  a  solid  right-circular  cylindrical  steel  projectile.  The  results  obtained  show  that 
irreversible  densification  can  provide  only  a  minor  improvement  in  the  ballistic  resistance  of 
glass  and  only  in  the  case  of  high-velocity  (ca.  1000  m/s)  projectiles.  Furthermore,  it  was  dem¬ 
onstrated  that  if  the  key  irreversible  compaction  parameters  can  be  adjusted  by  modifications  in 
glass  chemistry  and  microstructure,  significant  improvements  in  the  glass  ballistic  resistance  can 
be  attained  over  a  relatively  wide  range  of  projectile  velocities. 

Keywords:  glass,  molecular-level  modelling  and  simulations,  high-pressure  irreversible 
densification 


1  INTRODUCTION 

A  public  domain  literature  review  carried  out  as  part 
of  this  study  revealed  that  several  different  materials 
and  design  strategies  are  currently  being  used  in 
transparent  ballistic-impact  resistant  vehicle  struc¬ 
tures  (e.g.  windshields,  door  windows,  viewports, 
etc.).  Among  the  most  recently  introduced  transpar¬ 
ent  materials  and  technologies,  the  following  have 
received  the  most  attention:  transparent  crystalline 
ceramics  (e.g.  aluminium-oxynitride  spinel,  AlON, 
sapphire  [1]),  new  transparent  polymer  materials 
(e.g.  transparent  nylon  [2]),  new  interlayer  technolo¬ 
gies  (e.g.  polyurethane  bonding  layers  [3]),  and  new 
laminate  structure  designs  [4].  Despite  the  clear 
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benefits  offered  by  these  materials  and  technologies 
(e.g.  transparent  ceramics  offer  a  very  attractive  com¬ 
bination  of  high  stiffness  and  high  hardness  levels, 
highly  ductile  transparent  polymers  provide  superior 
fragment  containing  capabilities,  etc.),  ballistic  glass 
remains  an  important  constituent  material  in  a 
majority  of  transparent  impact  resistant  structures 
used  today.  Among  the  main  reasons  for  the  wide- 
scale  use  of  glass,  the  following  three  are  most  fre¬ 
quently  cited: 

(a)  glass-structure  fabrication  technologies  enable 
the  production  of  curved,  large  surface- area, 
transparent  structures  with  thickness  approach¬ 
ing  several  inches; 

(b)  relatively  low  material  and  manufacturing  costs; 

(c)  compositional  modifications,  chemical  strength¬ 
ening,  and  controlled  crystallization  have  dem¬ 
onstrated  to  be  capable  of  significantly 
improving  the  ballistic  properties  of  glass  [2] . 
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The  development  of  new  glass-based  transparent 
impact  resistant  structures  aimed  at  reducing  the 
vulnerability  of  protected  vehicle  occupants,  and 
on-board  instrumentation  to  various  threats  typi¬ 
cally  includes  extensive  prototyping  and  labora- 
tory/field  testing.  These  prototyping/testing 
programs  are  critical  for  ensuring  the  utility  and 
effectiveness  of  the  transparent  impact  resistant 
structures.  However,  the  use  of  prototyping/testing 
programs  is  generally  expensive,  time-consuming, 
and  involving  destructive  test  procedures.  While 
the  role  of  prototyping/testing  programs  remains 
critical,  they  are  increasingly  being  complemented 
by  the  corresponding  computation-based  modelling 
and  simulation  efforts.  However,  the  availability  of 
realistic  physically  based  material  models  describing 
deformation /fracture  response  of  ballistic  glass 
under  high-deformation-rate/high-pressure  loading 
conditions  is  one  of  the  key  requirements  for 
attaining  a  high  level  of  utility  and  fidelity  of 
these  computation-based  modelling  and  simulation 
approaches.  Therefore,  one  of  the  main  objectives  of 
this  study  is  to  further  advance  the  application 
of  computational  modelling/ simulation-based  engi¬ 
neering  approaches  of  transparent  impact-resistant 
structures  via  improvements  in  the  accuracy  of  the 
existing  ballistic-glass  material  models. 

A  comprehensive  literature  review  carried  out  as 
part  of  this  study  revealed  that  the  mechanical  behav¬ 
iour  of  glass  is  modelled  predominantly  using  three 
distinct  approaches: 

(a)  molecular-modelling  methods; 

(b)  continuum-material  approximations; 

(c)  models  based  on  explicit  crack  representation. 

A  brief  overview  and  the  main  findings  for  each  of 
these  three  approaches  are  given  in  the  remainder 
of  this  section. 

1.1  Molecular-level  material  modelling 

The  first  molecular-level  computational  investigation 
of  glass  reported  in  open  literature  can  be  traced 
back  to  the  1976  work  of  Woodcock  et  al.  [5],  Since 
that  time,  major  advances  in  computer  technol¬ 
ogy  and  the  introduction  of  high-fidelity  quantum 
mechanics-based  force  fields  (inter-atomic  poten¬ 
tials)  have  allowed  for  more  accurate  computational 
modelling  of  glass  elastic  constants,  strength,  chem¬ 
ical  and  thermal  diffusivities,  surface  energies,  etc. 
Of  interest  to  this  study,  a  great  number  of  research¬ 
ers  have  investigated,  using  molecular  modelling 
techniques,  the  propensity  of  various  types  of  glass 
(of  different  chemistries  and  microstructures) 
to  undergo  irreversible  (permanent)  densification 


when  subjected  to  high  hydrostatic  pressures  on  the 
order  of  10  GPa  [6-8].  The  emphasis  in  these  investi¬ 
gations  was  placed  on  elucidating  the  main  atomic- 
level  mechanisms  and  processes  (e.g.  increased 
coordination  number,  often  referred  to  as  coordina¬ 
tion  defects,  creation  of  new  metastable  chemical 
bonds,  etc.)  associated  with  high-pressure  irrevers¬ 
ible  densification  of  glass.  In  this  study,  on  the  other 
hand,  molecular-modelling  investigations  of  high- 
pressure  irreversible  densification  of  glass  will  be 
carried  out  in  order  to  assess  its  effect  on  the  conti¬ 
nuum-level  pressure  versus  degree-of-compression 
(the  negative  of  volumetric  strain)  relation,  also 
known  as  the  equation  of  state  (EOS).  In  addition, 
the  ability  of  glass  densification  to  act  as  a  potent 
energy  absorbing  process,  and  thus,  as  a  glass-tough¬ 
ening  mechanism  is  investigated. 

1.2  Continuum-level  material  modelling 

Within  the  continuum-level  glass  models  [9-15], 
glass  is  treated  as  a  continuum  material  whose  stiff¬ 
ness  and  strength  properties  may  become  degraded 
by  nucleation,  growth,  and  coalescence  of  cracks.  The 
fundamental  assumption  in  these  models  is  that  the 
elastic  stiffness  and  strength  degradations  are  the 
results  of  inelastic  deformation  caused  by  micron 
and  submicron  size  cracks,  and  that  this  degradation 
can  be  quantified  using  a  so-called  ‘damage  tensor’ 
whose  evolution  during  loading  can  be  formulated 
using  generalized  Griffith-type  crack  initiation 
and  propagation  criteria  for  brittle  materials.  In  addi¬ 
tion,  some  continuum  models  account  for  the 
interactions  between  the  cracks,  their  coalescence, 
friction  between  fragments,  competition  between 
micro -cracking  leading  to  fine-scale  fragmentation 
of  glass  and  macro -cracldng  giving  rise  to  coarse  frag¬ 
mentation,  etc.  In  addition  to  the  physically  based 
continuum-material  models  for  glass  mentioned 
above,  the  Johnson-Holmquist  2  (commonly  referred 
to  as  the  JH2  model  [16])  is  often  used  to  model  the 
behaviour  of  glass  under  high-loading  rate  condi¬ 
tions.  Despite  its  phenomenological  nature,  the 
JH2  model  has  been  found  to  often  provide  a  reason¬ 
ably  good  account  for  glass  response  under 
these  loading  conditions.  For  this  reason,  the  JH2 
model  will  be  used  in  this  study  for  coupling  with 
the  molecular  modelling  approaches  mentioned 
above.  Specifically,  the  polynomial  EOS  used  to 
account  for  the  hydrostatic /volumetric  response  of 
glass  within  the  JH2  model  will  be  modified  to 
include  pressure  versus  degree-of-compression 
results  which  will  be  obtained  using  the  aforemen¬ 
tioned  molecular  modelling  procedure.  In  addition, 
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molecular-level  results  pertaining  to  irreversible  den- 
sification-induced  strengthening  will  be  used  to 
modify  the  JH2  strength  constitutive  law. 

1.3  Explicit  crack  representation 
material  models 

Within  this  material-modelling  framework,  glass  is 
treated  as  a  linear  elastic  material,  and  its  fracture  is 
considered  to  take  place  via  nucleation,  propagation, 
and  coalescence  of  discrete  (rather  than  smeared- 
out/homogenized)  cracks  during  impact  [17].  In 
other  words,  while  within  the  continuum  modelling 
framework,  the  stiffness/strength-degrading  effect  of 
smeared-out  cracks  is  included  only  implicitly;  in  the 
explicit  crack  representation  material  models,  cracks 
are  considered  as  discrete  entities  and  their  effect  on 
material  stiffness /strength  is  accounted  for  explicitly. 
When  the  latter  type  of  glass  models  are  implemented 
into  a  finite- element  computational  framework, 
crack  nucleation  and  propagation  are  handled  by 
duplicating  nodes  at  the  crack  tip/front.  Adaptive 
re-meshing  is  used  to  provide  a  rich  enough  set  of 
possible  fracture  paths  around  the  crack  tip.  As  a 
crack  grows,  forces  at  newly  cracked  (free)  surfaces 
are  brought  to  zero  in  accordance  with  the  Griffith 
criterion  to  account  for  crack  growth  induced  unload¬ 
ing.  This  enables  explicit  modelling  of  the  crack 
coalescence  process  which  can  lead  to  fragment  for¬ 
mation.  The  major  disadvantages  of  the  discrete 
models  are  that  they  are  extremely  computationally 
expensive  and  become  intractable  as  the  number  of 
cracks  increases.  That  is,  in  order  to  capture  all  pos¬ 
sible  crack  nucleating  sites,  meshes  with  micron-size 
element  are  ultimately  required.  Hence,  despite  the 
fact  that  inclusion  of  high-pressure  densification 
effects  may  also  be  beneficial  to  the  discrete  glass 
models,  due  to  their  prohibitively  high  computational 
costs,  they  will  not  be  considered  in  this  study. 

As  demonstrated  above,  molecular-level,  conti¬ 
nuum-level,  and  discrete  modelling  are  maturing 
areas  of  glass  research  which  are  capable  of  revealing 
complex  intrinsic  mechanisms  and  phenomena  asso¬ 
ciated  with  deformation  and  fracture  in  glass. 
However,  these  modelling  approaches  are  typically 
concerned  only  with  the  effect  of  the  observed  pro¬ 
cesses/mechanisms  on  the  behaviour  of  glass  at  their 
respective  length  scales  and  practically  no  reports 
were  found  where  the  knowledge  about  glass  behav¬ 
iour  at  one  length  scale  was  used  to  improve  glass 
models  at  other  length  scale  (s).  Therefore,  the  main 
objectives  of  this  study  are: 

(a)  to  investigate  and  quantify  (using  molecular-level 
modelling  and  simulation  techniques)  the 


irreversible-densification  process  in  soda-lime 
glass  at  high  pressures; 

(b)  to  determine  if  modifications  in  the  EOS  and  the 
strength  constitutive  law  of  the  JH2  continuum- 
material  model  for  soda-lime  glass  to  include  the 
effects  of  high-pressure  irreversible  densification 
obtained  in  (a)  significantly  alters  the  mechanical 
response  of  glass  under  ballistic-loading 
conditions; 

(c)  to  carry  out  a  preliminary  assessment  of  high- 
pressure  irreversible  densification  as  an  energy¬ 
absorbing/ strength-enhancing  mechanism  in 
glass. 

Towards  that  end,  molecular- level  calculations  are 
carried  out  first  to  quantify  the  basic  pressure 
versus  degree-of-compression  relation  in  glass  at 
high-pressure  as  this  material  undergoes  irreversible 
densification.  Also,  molecular-level  simple  shear  tests 
were  carried  out  to  assess  the  extent  of  irreversible 
densification-induced  strengthening.  The  results 
obtained  are  next  used  to  modify  the  JH2  EOS  and 
strength  constitutive  law  for  glass.  Then,  a  series  of 
transient  non-linear  dynamics  analyses  of  transverse 
impact  of  a  glass  test  panel  with  a  solid  right-circular 
cylinder  fragment  simulating  projectile  (FSP)  are  car¬ 
ried  out  to  examine  (a)  the  extent  of  change  in  the 
temporal  and  spatial  distribution  of  deformation 
and  damage  and  (b)  the  resulting  change  in  ballis¬ 
tic-penetration  resistance  brought  about  by  the 
aforementioned  modifications  in  the  glass  material 
model. 

The  organization  of  this  article  is  as  follows:  a  dis¬ 
cussion  of  the  inter-atomic  force  field  potentials, 
computational  cell,  and  the  computational  method 
used  in  the  molecular-level  simulations  are  all  pre¬ 
sented  in  section  2.  Then,  a  brief  overview  of  the 
JH2  continuum-material  model  for  glass  is  provided 
in  section  3.  The  main  results  obtained  in  this  study 
(including  the  definition  of  a  FSP  transverse-impact 
problem  for  a  glass  test-panel)  are  presented  and  dis¬ 
cussed  in  section  4,  while  the  key  conclusions  result¬ 
ing  from  this  study  are  summarized  in  section  5. 

2  MOLECULAR-LEVEL  MODELLING  OF  GLASS 

At  the  molecular  level,  soda-lime  glass  is  treated  as  a 
discrete  material  consisting  of: 

(a)  silicon  (Si)  and  oxygen  (O)  atoms  mutually 
bonded  via  a  single  covalent  bond  and  forming 
a  connected,  non-structured/amorphous  net¬ 
work  of  silica  (SiO^)  tetrahedra; 

(b)  oxygen  anions  (O2-)  attached  as  terminal  func¬ 
tional-groups  to  the  fragmented  silica  tetrahedra 
network; 
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(c)  sodium  cations  (Na+)  dispersed  between  frag¬ 
mented  silica  tetrahedra  networks  and  ionically 
bonded  to  the  oxygen  anions. 

To  fully  account  for  the  bonding  and  non-bonding 
types  of  interactions  between  the  atoms /ions  listed 
above,  one  must  define  the  respective  interaction- 
potential  functions  (commonly  referred  to  as  the 
force-fields),  as  well  as  the  associated  atomic-polar 
and  ionic  charges. 

While  glass  is  an  amorphous  material  and  does  not 
possess  any  long-range  regularity  in  its  atomic/ 
molecular  structure,  modelling  of  bulk  behaviour  of 
glass  is  typically  done  at  the  molecular  level  by 
assuming  the  existence  of  a  larger  unit  cell. 
Repetition  of  this  cell  in  the  three  orthogonal  direc¬ 
tions  (the  process  also  known  as  application  of  the 
‘periodic  boundary  conditions’)  results  in  the  forma¬ 
tion  of  an  infinitely  large  bulk-type  material. 

Molecular-modelling  simulations  typically  rely  on 
one  of  the  following  two  techniques: 

(a)  molecular  statics,  the  technique  within  which  the 
potential  energy  of  the  molecular  structure  in 
question  is  minimized  with  respect  to  the  posi¬ 
tion  of  the  constituent  and  ions  in  the  unit  cell  as 
well  as  with  respect  to  the  size  and  shape  of  the 
unit  cell; 

(b)  molecular  dynamics,  a  technique  within  which 
the  Newton’s  equations  of  motion  are  solved  for 
all  interacting  atoms  an  ions  in  the  system  as  a 
function  of  time  and  the  appropriate  ensemble 
averages  are  used  to  assess/quantify  various 
molecular-level  material  properties. 

The  discussion  presented  above  indicates  that  the 
three  main  components  of  a  molecular-level  model 
which  must  be  defined  are: 

(a)  the  force-fields; 

(b)  the  initial  unit  cell  size  and  shape  as  well  as  the 
initial  positions  of  atoms  and  ions  within  it; 

(c)  the  computational  procedure  to  be  employed. 

A  brief  description  of  these  three  components  used  in 
the  present  molecular- level  investigation  of  soda- 
lime  glass  behaviour  under  high  pressures  is  pre¬ 
sented  in  the  remainder  of  this  section. 


2.1  Force-fields 

While  accurate  simulations  of  a  system  of  interacting 
particles  (i.e.  atoms  or  ions)  generally  entail  the  appli¬ 
cation  of  quantum  mechanical  techniques,  such 
techniques  are  computationally  quite  expensive  and 
are  usually  feasible  only  in  systems  containing  up  to  a 
few  hundreds  of  interacting  particles.  In  addition,  the 


main  goal  of  simulations  of  the  systems  containing  a 
large  number  of  particles  is  generally  to  obtain  the 
systems’  bulk  properties,  which  are  primarily  con¬ 
trolled  by  the  location  of  atomic  nuclei  and  the 
knowledge  of  the  electronic  structure,  provided  by 
the  quantum  mechanics  techniques,  is  not  critical. 
Under  these  circumstances,  a  good  insight  into  the 
behaviour  of  a  system  can  be  obtained  if  a  reasonable, 
physically  based  approximation  of  the  potential 
(force-field)  in  which  atomic  nuclei  move  is  available. 
Such  a  force-field  can  be  used  to  generate  a  set  of 
system  configurations  which  are  statistically  consis¬ 
tent  with  a  fully  quantum-mechanical  description. 

As  stated  above,  a  crucial  point  in  the  molecular- 
level  simulations  of  multi-particle  systems  is  the 
choice  of  the  force-fields  which  describe,  in  an 
approximate  manner,  the  potential  energy  hyper-sur¬ 
face  on  which  the  atomic  nuclei  move.  In  other 
words,  the  knowledge  of  force-fields  enables  determi¬ 
nation  of  the  potential  energy  of  a  system  in  a  given 
configuration.  In  general,  the  potential  energy  of  a 
system  of  interacting  particles  can  be  expressed  as  a 
sum  of  the  valence  (or  bond),  Ev;ih.nc(.}  cross-term, 
^cross-term,  and  non-bond,  EnonAmn(h  interaction  ener¬ 
gies  as 

Etotal  =  ^valence  "T  Across— term  T  ^non-bond  ( 1 ) 

The  valence  energy  generally  includes  a  bond  stretch¬ 
ing  term,  £:bond,  a  two-bond  angle  term,  Eangle,  a  dihe¬ 
dral  bond-torsion  term,  Etorsion,  an  inversion  (or  an 
out-of-plane  interaction)  term,  Eoop,  and  a  Urey- 
Bradlay  term  (which  involves  interactions  between 
two  particles  bonded  to  a  common  particle),  f:UB,  as 

Eyalence  —  Ebond  T  Wangle  T  -Etorsion  ~t"  Eoop  +  Eub  (2) 

A  schematic  explanation  of  the  first  four  types  of 
valence  atomic  interactions  is  given  in  Fig.  1 . 

The  cross-term  interacting  energy,  Ecross_term, 
accounts  for  the  effects  such  as  bond  length  and 


(a)  ,  (b) 

b 


Fig.  1  A  schematic  of  (a)  stretch,  (b)  angle,  (c)  torsion, 
and  (d)  inversion  valence  atomic  interactions 
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angle  changes  caused  by  the  surrounding  atoms  and 
generally  includes: 

(a)  stretch-stretch  interactions  between  two  adja¬ 
cent  bonds,  £bond— bond* 

(b)  stretch-bend  interactions  between  a  two-bond 
angle  and  one  of  its  bonds,  £bond-angie; 

(c)  bend-bend  interactions  between  two  valence 
angles  associated  with  a  common  vertex  particle, 

E^angl  e-angle, 

(d)  stretch-torsion  interactions  between  a  dihedral 
angle  and  one  of  its  end  bonds,  Tendjbond-torsionl 

(e)  stretch-torsion  interactions  between  a  dihedral 
angle  and  its  middle  bond,  -Emiddie_bond-torsion> 

(f)  bend-torsion  interactions  between  a  dihedral 
angle  and  one  of  its  valence  angles,  Ekngie-torsion; 

(g)  bend-bend-torsion  interactions  between  a 
dihedral  angle  and  its  two  valence  angles; 

Eiangl  e-angle-torsion, 

terms  as 


-Ecross— term  —  Ebond— bond  T  EaIlg|e_nng]e  +  Ebond— angle 
~b  E’endjjond-torsion  T  Ebii(j(iie_b0ncl — torsion 
"b  E'angie_torsion  +  £angle— angle— torsion 

(3) 


The  non-bond  interaction  term,  Enon-hond,  accounts 
for  the  interactions  between  non-bonded  particles 
and  includes  the  van  der  Waals  energy,  Ev,lw,  the 
Coulomb  electrostatic  energy,  ficouiomb,  and  the 
hydrogen  bond  energy,  £H_bond.  as 

Enon— bond  —  EvdVV  ~b  Edoulomb  "b  Eh— bond  (4) 


Particle  interactions  in  the  soda-lime  glass  system 
under  investigation  are  modelled  using  COMPASS 
(condensed-phased  optimized  molecular  potential 
for  atomistic  simulation  studies),  the  first  ah  initio 
force-field  that  enables  an  accurate  and  simultaneous 
prediction  of  various  gas-phase  and  condensed- 
phase  properties  of  organic  and  inorganic  materials 
[18,  19],  The  COMPASS  force-field  uses  the  following 
expression  for  various  components  of  the  potential 
energy 

Ebond  =  X  [K^b  -  b0)2+K3(b  ~  hf+^b  -  b0 )4] 
b 

(5) 

Eangle  =  X  [H2@  ~  6o)2+H3(9  -  Bnf+HA(9  -  0O)4] 

e 

(6) 

Etorsion  =  ^  '  [Vl[l  —  COS  (0  — 

<P 

+  V^fl  —  cos  (20  —  0^)] 

+  V3[l-  cos  (30-0°)]]  (7) 


Eoop  =  X 

a: 

(8) 

-^bond— bond  —  ^  ^  M) ){b  ~  ^o) 

(9) 

b  b' 

^angle-angle  =  X^  XI  Fq&(0  ~  #0 )(&  ~  O'q) 

e  o' 

(10) 

-^bond— angle  ~  ^  ^  ]  Ffo q{1)  —  bo)(0  —  Oq) 

b  e 

(11) 

^end_bond— torsion  —  ^  ^  ]  Fb(f){b  —  bo)[Vi  COS  0 

b  ip 


+  V2  cos  20  +  V3  cos  30]  (12) 

-^middle.bond— torsion  — 

y  <p 

x  [Fi  cos  0  +  F2  cos  20  +  F3  cos  30] 

(13) 


Eangle— torsion  —  EE  Fo(j)(6  —  Oq) 
e  0 

[Vi  cos  0  +  V2  cos  20+  V3cos30]  (14) 

-Wangle— angle— torsion  —  EEE  Kpee' 

(pee1 

x  cos  0(0  —  E0)(ff  —  $(,)  (15) 

(,i(,i  (16) 


■^Coulomb 


_  HiH] 
~  “  era 


i>j 


F  _  \ "  Aii  B’i 
bvdW  —  2_,  ,-9  „6 

i>j  \_H  ij_ 


(17) 


where  b  and  b  are  the  bond  lengths,  9  the  two-bond 
angle,  <p  the  dihedral  torsion  angle,  /  the  out  of  plane 
angle,  q  the  atomic  charge,  e  the  dielectric  constant, 
and  r,j  Lhe  i-j atomic  separation  distance.  b0>  Kt  (i—2- 
4),  0O,  Hi  ( j  =  2-4),  $  (i=  1-3),  Vt  (/=  1-3),  Fw,  b'n,  Fm, 
0q  ,  E/x,,  Fhrp  Ffyqj,  F[  ( /  —  1—3),  Ffj(p>  KyjQQ* ,  Ay,  and  By  are 
the  material-system  dependent  parameters  imple¬ 
mented  into  Discover  [20],  the  atomic  simulation 
program  used  in  this  study. 


2.2  Computational  cell 


In  the  molecular-level  modelling  portion  of  this 
study,  a  2916-particle  computational  cell  was  created 
with  a  chemical  composition  of  (Na2O)0.i5  (SiO2)0.85- 
The  molecular-level  unit  cell  was  constructed  with  a 
cubic  geometry  having  a  uniform  edge-length  of 
3.347  nm.  An  amorphous  cell  of  sodium  silicate 
glass  was  thus  created  with  a  density  of  2.613  g/cm3. 
The  three  edges  (a,  b,  and  c)  of  the  cell  were  aligned, 
respectively,  with  the  three  coordinate  axes  (x,  y,  and 
z).  To  create  a  bulk-material  environment  for  the  par¬ 
ticles,  periodic  boundary  conditions  were  applied  to 
all  six  faces  of  the  unit  cell. 
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Fig.  2  The  computational  unit  cell  for  soda-lime  glass 
molecular-level  simulations  used  in  this  study 


To  create  the  initial  particle  configuration  in  the 
unit  cell,  the  Visualizer  [21]  from  Accelrys  was  first 
used  to  construct  a  short  silica-chain  fragment.  The 
fragment  was  then  ‘grown’  by  a  duplicate-and-attach 
process  using  the  same  program.  The  resulting  silica 
network  was  next  used  within  the  Amorphous  Cell 
program  [22]  from  Accelrys  to  randomly  populate 
the  computational  cell  while  ensuring  that  the 
target  material  density  of  2.613  g/cm3  was  attained. 
Finally,  sodium  cations  were  added  to  obtain  an  elec¬ 
trically  neutral  system  of  particles.  An  example  of  a 
typical  molecular- level  topology  within  a  single  unit 
cell  is  displayed  in  Fig.  2.  When  creating  the  compu¬ 
tational  cell  for  the  subject  material  (soda-lime  glass 
in  the  present  case),  one  should  try  to  answer  the 
question  as  to  how  representative  is  the  computa¬ 
tional  cell  of  the  material  in  question.  Typically,  the 
size  of  the  computational  cell,  i.e.  the  number  of  the 
constituent  atoms  analysed  is  a  compromise  between 
the  available  computer  resources  (CPU  time,  primar¬ 
ily)  and  a  desire  to  maximize  the  size  of  the  compu¬ 
tational  model.  The  reason  that  the  computational- 
model  size  should  be  maximized  is  that  the  effect 
of  the  periodic  boundary  conditions  which  artificially 
introduce  configurational-order  into  an  amorphous/ 
disordered  solid  should  be  minimized.  To  qualify  the 
computational  cell  used  in  this  study,  few  simulations 
are  carried  out  using  a  computational  cell  with  a 
doubled-edge  length.  Since  no  statistically  significant 
differences  in  the  key  results  were  seen  for  the  two 
choices  of  the  unit  cell,  it  is  concluded  that  the  com¬ 
putational  cell  used  is  adequate. 

2.3  Computational  method 

High-pressure  irreversible  densification  of  glass  was 
studied  using  a  conventional  NPT  dynamics  method 
within  the  Discover  program  [20]  from  Accelrys, 


where  N  (the  number  of  particles),  P  (pressure),  and 
T  (temperature)  are  the  system  variables  that  are  held 
constant  or  ramped  in  a  controlled  manner  during 
compression-simulation  runs.  The  equations  of 
motion  were  integrated  using  the  velocity  Verlet  algo¬ 
rithm  with  a  time  step  of  1 .0  fs.  To  comply  with  typical 
ballistic-impact  loading  conditions,  pressure  was 
ramped  linearly  at  a  rate  of  0.3GPa/ps.  Pressure 
ramping  was  accomplished  by  3  GPa  discrete  incre¬ 
ments  in  pressure  followed  by  10.0  ps  equilibration 
times  at  a  given  pressure  level.  Temperature,  on  the 
other  hand  was  held  constant  at  a  value  of  300  K  using 
a  Nose  Thermostat  [23]. 

At  the  end  of  equilibration  at  each  pressure  level, 
the  average  material  density  was  computed  from  the 
corresponding  pair- correlation  functions.  This  pro¬ 
cedure  yielded  the  sought  pressure  versus  degree- 
of-compression  relation  for  the  pressurization  por¬ 
tion  of  a  loading  cycle.  To  determine  the  correspond¬ 
ing  pressure  versus  degree-of-compression  relation 
during  the  depressurization  portion  of  a  loading 
cycle,  pressure  was  decreased  in  a  similar  manner 
as  discussed  above.  Differences  in  the  pressure 
versus  degree-of-compression  relations  for  the  pres¬ 
surization  and  depressurization  portions  of  the 
loading  cycle  are  then  used  to  quantify  the  extent  of 
high-pressure  irreversible  densification  of  glass. 

As  will  be  explained  later  in  more  detail,  irreversible 
densification  can  affect  the  ballistic-penetration  of 
glass  not  only  via  increased  density  but  also  through 
the  associated  changes  in  the  material  strength.  To 
assess  the  extent  of  densification-induced  strength¬ 
ening,  molecular-level  simple-shear  tests  were  car¬ 
ried  out.  These  tests  were  conducted  through  the 
use  of  a  Discover  input  file  which  was  written  in  a 
basic  tool  command  language  (BTCL).  This  enabled 
the  use  of  a  scripting  engine  that  provides  very  precise 
control  of  simulation  runs,  e.g.  a  cell  deformation  to 
be  carried  out  in  small  steps  each  followed  by  energy 
minimization.  The  minimization  portion  of  the  mole¬ 
cular-level  simple  shear  tests  was  carried  out  using  a 
combination  of  three  (steepest  descent,  conjugate 
gradient,  and  Newton’s)  potential- energy  minimiza¬ 
tion  algorithms  within  Discover  [20].  These  algo¬ 
rithms  are  automatically  activated/deactivated  as 
the  molecular-level  configuration  approaches  its 
energy  minimum  (i.e.  the  steepest  descent  method 
is  activated  at  the  beginning  of  the  energy- minimiza¬ 
tion  procedure,  while  the  Newton’s  method  is  utilized 
in  the  past  stages  of  the  simulation). 

3  JH2  CONTINUUM-MATERIAL  MODEL 

To  completely  define  a  continuum-level  material 
dynamic  model,  the  relationships  between  the  flow 
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variables  (pressure,  mass-density,  energy- density, 
temperature,  etc.)  must  be  specified.  These  relations 
typically  involve: 

(a)  an  EOS; 

(b)  a  strength  equation; 

(c)  a  failure  equation; 

(d)  an  erosion  equation. 

These  equations  arise  from  the  fact  that,  in  general, 
the  total  stress  tensor  can  be  decomposed  into  a 
sum  of  a  hydrostatic  stress  (pressure)  tensor 
(which  causes  a  change  in  the  volume /density  of 
the  material)  and  a  deviatoric  stress  tensor  (which 
is  responsible  for  the  shape  change  of  the  material) . 
An  EOS  then  is  used  to  define  the  corresponding 
functional  relationship  between  pressure,  mass  den¬ 
sity  (degree-of-compression),  and  internal  energy 
density  (temperature).  Likewise,  a  strength  relation 
is  used  to  define  the  appropriate  equivalent  plastic 
strain,  equivalent  plastic  strain  rate,  and  tempera¬ 
ture  dependencies  of  the  materials  yield  strength. 
This  relation,  in  conjunction  with  the  appropriate 
yield- criterion  and  flow-rule  relations,  is  used  to 
compute  the  deviatoric  part  of  stress  under  elastic- 
plastic  loading  conditions.  In  addition,  a  material 
model  generally  includes  a  failure  criterion  (i.e.  an 
equation  describing  the  hydrostatic  or  deviatoric 
stress  and/or  strain  condition(s)  which,  when  met, 
causes  the  material  to  fracture  and  lose  its  ability  to 
support,  abruptly  in  the  case  of  brittle  materials  or 
gradually  in  the  case  of  ductile  materials,  normal  and 
shear  stresses).  Such  a  failure  criterion  in  combination 
with  the  corresponding  material-property  degrada¬ 
tion  and  the  flow- rule  relations  governs  the  evolution 
of  stress  during  failure.  The  erosion  equation  is  gen¬ 
erally  intended  for  eliminating  numerical  difficulties 
arising  from  highly  distorted  Lagrange  cells. 
Nevertheless,  the  erosion  equation  is  often  used  to 
provide  additional  material  failure  mechanisms, 
especially  in  materials  with  limited  ductility. 

To  summarize,  the  EOS  along  with  the  strength  and 
failure  equations  (as  well  as  with  the  equations  gov¬ 
erning  the  onset  of  plastic  deformation  and  failure 
and  the  plasticity  and  failure  induced  material  flow) 
enable  assessment  of  the  evolution  of  the  complete 
stress  tensor  during  a  transient  non-linear  dynamics 
analysis.  Such  an  assessment  is  needed  in  order  to 
solve  the  governing  (mass,  momentum,  and  energy) 
conservation  equations.  It  is  important  to  note  that 
separate  evaluations  of  the  pressure  and  the  deviato¬ 
ric  stress  enable  inclusion  of  the  non-linear  effects  in 
the  EOS.  Generally,  these  effects  are  shock  related 
but,  in  this  study,  they  will  be,  at  least  partly,  attrib¬ 
uted  to  the  phenomenon  of  high-pressure  irreversible 
densification. 


In  this  study,  glass  was  modelled  using  the  JH2 
brittle-material  model  [24,  25].  The  JH2  model  is  a 
phenomenological  model  which  postulates  the  exis¬ 
tence  of  two  terminal  glass  states:  (a)  an  intact  mate¬ 
rial  and  (b)  a  failed  material.  The  two  material  states 
are  weighted  by  a  single  scalar  variable  called 
damage,  D,  whose  evolution  is  governed  by  an  inelas¬ 
tic  (plasticity-like  deformation  model).  The  JH2 
model  includes  a  polynomial-type  EOS,  a  strength 
model  (based  on  the  von  Mises  yield  criterion,  nor¬ 
mality  flow  rule  and  a  pressure  and  strain-rate  hard¬ 
ening  constitutive  relation),  a  progressive  failure 
model,  and  an  instantaneous  geometric  strain- 
based  erosion  criterion.  The  values  of  all  the  JH2 
material-model  parameters  for  soda-lime  float  glass 
are  available  in  the  ANSYS/Autodyn  materials  library 
[26].  Further  details  of  the  JH2  model  for  brittle  mate¬ 
rials  are  provided  in  the  remainder  of  this  section. 

3.1  Polynomial  EOS 

Within  the  JH2  polynomial  EOS,  the  effect  of  internal- 
energy  density  is  neglected  and  the  pressure  versus 
degree-of-compression  in  a  damage-free  material  is 
defined  as 

P  =  Kiii  -T  Kip.2  4"  Ki/j.3 ,  /1  >  0 

(hydrostatic  compression)  (18) 

and 

P  =  K\ji ,  fi  <  0  (hydrostatic  tension)  (19) 

where  degree-of-compression  is  p,  =  (p/p0  -  1)  and  p 
the  current  density,  while  p0  (the  reference  density), 
Ki  (the  bulk  modulus),  and  K2  and  K3  the  material- 
specific  constants. 

After  glass  has  begun  to  accumulate  damage  (i.e. 
when  the  extent  of  damage  is  no  longer  zero,  D  >  0), 
equation  (18)  has  to  be  upgraded  to  include  the  effect 
of  bulking.  Bulking  is  a  phenomenon  associated  with 
the  fact  that  fragments  of  fractured  materials  are  not 
generally  fully  conformable  and,  consequently  frac¬ 
tured  material  is  associated  with  a  larger  volume  (a 
lower  density  at  a  constant  pressure)  than  the 
damage-free  material.  The  bulking  modified  polyno¬ 
mial  EOS  is  then  given  by  [26] 

P  —  Kip  +  K2112  +  K^/i 3  +  A P,  /i  >  0  (20) 

where  the  hulking-induced  pressure  increment,  A P, 
is  determined  from  energy  considerations  and  varies 
from  zero  at  D  =  0  to  A  Pmax  at  D  =  1 .0.  Assuming  that 
a  fraction  of  the  internal  elastic  energy  decrease  (due 
to  decrease  in  deviatoric  stresses  in  the  material)  is 
converted  to  an  increase  in  potential  internal  energy, 
the  bulking  induced  pressure  increment  A  Pat  a  time 
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t  +  At  can  be  represented  in  terms  of  A P  at  the 
time  fas 


AP(t  +  At)  —  -Kfi(t  +  At) 


I  (Ki/j,(t  +  At)  +  AP(t))2 
+2fiK\  A  U) 


(21) 


where  AU  is  the  decrease  in  deviatoric  elastic 
energy  due  to  damage  induced  yield- strength  reduc¬ 
tion  and  jJ>  the  fraction  of  the  deviatoric  elastic 
energy  converted  to  hydrostatic  potential/elastic 
energy.  The  decrease  in  deviatoric  elastic  energy  is 
given  by 

AU  —  Ut  —  Ut+At  (22) 

where 


Ut  = 


6  G 


(23) 


The  parameters  at  and  G  appearing  in  equation  (23) 
are  the  actual  yield  strength  and  the  damage-free 
shear  modulus  of  the  glass  material,  respectively. 


3.2  Strength  model 

Within  the  JH2  strength  model  [24,  25] ,  the  normal¬ 
ized  yield  strength  is  defined  as  a  damage -weighted 
rule  of  mixtures  of  the  corresponding  damage-free 
and  fractured  yield  strengths  as 

a*  —  a*  —  D(oj  -  a*f)  (24) 

where  subscripts  i  and /are  used  to  denote  intact  and 
fractured  material  states  and  the  superscript  *  indi¬ 
cates  that  the  corresponding  yield  strength  is  normal¬ 
ized  by  the  Hugoniot  elastic  limit  (HEL;  uni-axial 
stress)  yield  strength,  that  is 

a*  =  —  (25) 

cthel 

The  normalized  yield  strengths  cr*  and  aj  are  also 
defined  in  the  same  manner  as  a*.  The  normalized 
(pressure  and  strain  rate  dependent,  ideal-plastic) 
yield  strength  of  the  damage-free  material,  a*,  and 
the  fractured  material,  aj,  are,  respectively,  given  by 


a*  =  A(P*  +  T*)n(  1  +  Cine*) 


(26) 


Fig.  3  A  schematic  of  the  JH2  strength  model  in  the 
normalized  yield  strength,  cr*,  versus  normal¬ 
ized  pressure,  P  ,  plane 


and 


T 

■Phel 


(29) 


where  P  and  T  are  the  actual  pressure  and  the  maxi¬ 
mum  hydrostatic  tensile  pressure  that  the  glass  mate¬ 
rial  can  withstand,  respectively,  and  PHel  the  pressure 
at  HEL.  A  schematic  of  the  JH2  strength  model  in  the 
normalized  yield  strength  versus  normalized  pressure 
plane  is  displayed  in  Fig.  3. 

As  shown  in  reference  [25],  Phel  and  <tHel  are 
related  to  the  (uni-directional  shockwave-based)  uni- 
axial-strain  compressive  strength,  HEL,  as: 
HEL  =  PHel+  !°hel-  Since  both  PHel  and  <tHel  are 
dependent  on  the  compression  ratio  at  the  HEL, 
both  these  parameters  can  be  determined  from 
HEL.  Finally,  the  dimensionless  material  strain  rate, 
e*,  appearing  in  equations  (28)  and  (29)  is  defined  as 


(30) 


where  e  is  the  actual  strain  rate  and  e0  the  reference 
strain  rate  (typically  set  to  1.0  s-1). 


3.3  Failure  model 


Within  the  JH2  failure  model  [24,  25] ,  the  evolution/ 
accumulation  of  damage  is  defined  as 


ETp 


(31) 


of  =  B(P*)m(1  +  Cine*)  (27) 

where  A ,  B,  C,  M,  N,  and  ctHel  appearing  in  equations 
(25)  to  (27)  are  all  material- specific  constants  while  P 
and  T  are,  respectively,  defined  as 

P*=-^~  (28) 

Fuel 


where  Aep  is  the  increment  in  equivalent  plastic 
strain  with  an  increment  in  loading  and  the  failure 
strain  e£  is  a  pressure- dependent  equivalent  fracture 
strain,  which  is  defined  as 

efp  =  D\(P*  +  T*)m  (32) 

where  //  and  D2  are  material-specific  constants. 
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Within  the  JH2  failure  model,  fracture  occurs  when 
either  damage  reaches  a  critical  value  of  1.0  or  when 
negative  pressure  reaches  a  value  of  T.  Fracture 
material  has  no  ability  to  support  any  negative  pres¬ 
sure,  while  its  ability  to  support  shear  is  defined  by 
equation  (27). 

3.4  Erosion  model 

Within  the  finite-element  computational  frame¬ 
work,  numerical  difficulties  may  arise  from  exces¬ 
sive  distortion  of  the  elements.  To  overcome  these 
difficulties,  an  erosion  algorithm,  which  at  a  prede¬ 
fined  level  of  strain  removes  the  excessively  dis¬ 
torted  elements,  while  transferring  the  momentum 
associated  with  the  removed  nodes  to  the  remain¬ 
ing  nodes.  Following  our  prior  work  [27-30],  the 
erosion  criterion  is  defined  by  prescribing  a  critical 
value  for  the  instantaneous  geometrical  equivalent 
strain. 

4  RESULTS  AND  DISCUSSION 

4.1  Molecular-level  analysis  of  high-pressure 
densification  of  glass 

In  this  section,  a  brief  summary  of  the  molecular-level 
calculations  of  glass  response  to  high  pressure  is  pro¬ 
vided.  While  this  portion  of  the  work  yielded  numer¬ 
ous  results,  only  the  ones  directly  related  to  the 
potential  role  of  high-pressure  irreversible  densifica¬ 
tion  in  improving  impact/penetration  resistance  of 
glass  are  presented  and  discussed  in  greater  detail. 
Other  results  are  discussed  only  qualitatively. 

4.1.1  Molecular-level  topology 

A  detailed  examination  of  the  molecular-level  topol¬ 
ogy  after  subjecting  the  unit-cell  to  high  pressures 
revealed  distinct  differences  depending  on  whether 
the  maximum  pressure  was  below  or  above  ca.  4  GPa. 

Pressures  below  ~4GPa:  Molecular  modelling  of 
glass  pressurization/depressurization  revealed  that 
when  glass  is  exposed  to  pressures  not  exceeding  ca. 
4  GPa,  no  detectable  irreversible  changes  generally 
take  place  in  its  molecular  topology.  Closer  examina¬ 
tion  of  the  atomic  structure  at  different  pressure 
levels  between  0  and  4  GPa  revealed  as  follows. 

1.  The  presence  of  ‘active  regions’  within  which 
atoms  may  occasionally  undergo  large  displace¬ 
ments/jumps  (ca.  0.1  nm).  These  atomic  displace¬ 
ments  (the  results  not  shown  for  brevity)  were 
found  to  involve  coordinated  motion  of  at  least  a 
dozen  atoms  and  to  be  accompanied  by  abrupt 
changes  in  the  average  potential  energy. 


2.  In  most  cases,  atomic  rearrangement  described  in 

(a)  appears  to  be  associated  with  low-frequency 
transition  of  the  active  regions  between  two  dis¬ 
tinct  ‘equilibrium’  states  (of  comparable  potential 
energy) .  These  findings  are  in  complete  agreement 
with  those  found  by  Trachenko  and  Dove  [7]  who 
termed  this  phenomenon  as  Double  Well  Potential 
and  the  associated  low-frequency  transition/ 
vibrational  mode  as  the  ‘floppy  mode’. 

3.  While  pressurization  up  to  4  GPa  did  not  yield  any 
permanent  changes  in  the  molecular  topology,  the 
locations  of  the  active  regions  were  found  to 
change  with  pressure.  That  is,  the  regions  active 
at  one  pressure  level  may  become  inactive  at 
another  pressure  level  while,  at  the  same  time 
other  previously  inactive  regions  would  become 
active. 

Pressures  above  ~4  GPa:  When  the  computational 
cell  is  subjected  to  pressure  exceeding  ca.  4  GPa  and 
subsequently  depressurized  to  zero  pressure,  perma¬ 
nent  changes  in  the  glass  molecular  topology  were 
normally  observed.  These  molecular  topology 
changes  were  accompanied  by  a  permanent  density 
increase  in  the  order  of  3-7  per  cent.  An  example  of 
typical  results  obtained  in  this  portion  of  the  work  is 
given  in  Figs  4(a)  to  (d),  where  oxygen  atoms/anions 
are  displayed  in  red  and  silicon  atoms  are  shown  in 
grey  (as  well  green,  pink,  and  yellow  highlighting), 
while  sodium  cations  are  omitted  for  clarity.  To  aid 
in  visualization/interpretation  of  the  topological 
changes  experienced  by  glass  during  high-pressure 
loading/unloading  cycles,  only  a  30-40-atom  exem¬ 
plary  region  of  computational  cell  was  monitored  in 
Figs  4(a)  to  (d).  The  molecular  level  topologies  dis¬ 
played  in  these  figures  pertain  respectively  to: 

(a)  high-pressure  state,  Fig.  4(a); 

(b)  high-pressure  state  after  a  prolonged  (20  ps) 
relaxation  period,  Fig.  4(b); 

(c)  ambient  pressure  state  resulting  from  depressur¬ 
ization  of  state  (a),  Fig.  4(c); 

(d)  ambient  pressure  state  resulting  from  depressur¬ 
ization  of  state  (b),  Fig.  4(d). 

It  should  be  noted  that  the  results  displayed  in  Fig. 
4(c)  are  essentially  identical  to  those  in  the  initial  con¬ 
figuration  prior  to  loading.  A  closer  examination  of 
the  molecular- level  topology  results  displayed  in 
these  figures  revealed  as  follows. 

1.  Pressurization  alters  not  only  molecular-level 
topology,  but  also  changes  the  bonding  structure 
and  increases  the  average  coordination  number  (of 
mainly  Si  atoms).  This  can  be  seen  by  comparing 
the  results  displayed  in  Figs  4(a)  and  (c).  It  should 
be  noted  here,  that  the  results  displayed  in  Fig.  4(c) 
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Fig.  4  The  molecular  level  topologies  pertaining  to:  (a)  high-pressure  state;  (b)  high-pressure  state 
after  a  prolonged  (20  ps)  relaxation  period;  (c)  ambient  pressure  state  resulting  from  depres¬ 
surization  of  state  (a) ;  and  (d)  ambient  pressure  state  resulting  from  depressurization  of  state 
(b).  See  text  for  explanation 


are  used  in  place  of  the  initial  molecular-level 
topology  results.  In  these  figures,  it  is  seen  that  sil¬ 
icon  atoms  labelled  A  and  B  change  their  fourfold 
coordination  to  fivefold  coordination  upon 
pressurization. 

2.  As  implied  earlier,  if  the  depressurization  is  carried 
out  without  allowing  the  material  to  relax  at 
high  pressures,  the  molecular  level  configuration 
obtained  at  the  ambient  pressure  is  effectively 
identical  to  the  initial  configuration,  Fig.  4(c). 
Thus,  in  the  resultant  ambient-pressure  configura¬ 
tion,  most  silicon  atoms  regain  their  fourfold 
coordination. 

3.  Increased  duration  of  the  exposure  of  glass  to  high 
pressure  results  in  continued  changes  in  the  mole¬ 
cular-level  topology  and  bond  structure.  This  can 
be  seen  by  comparing  the  results  displayed  in 
Figs  4(a)  and  (b).  These  figures  show  that  relaxa¬ 
tion  of  glass  leads  to  the  C-label  silicon  atom 
acquiring  a  fivefold  coordination  while,  at  the 


same  time  smaller  size  Si-0  rings  are  being 
formed.  For  example,  a  twofold  ring  is  formed 
involving  the  A  and  C  silicon  atoms,  while  a  three¬ 
fold  ring  involving  the  C-,  D-,  and  ^-labelled  Si 
atoms  also  appears,  Fig.  4(b).  These  changes  in 
the  molecular-level  topology  of  glass  are  a  mani¬ 
festation  of  its  relaxation  to  a  material  state  that  is 
energetically  preferred  at  high  pressures. 

4.  Upon  depressurization  of  glass  which  was  relaxed 
at  high  pressures,  some  changes  in  the  molecular- 
level  topology  and  the  bonding  structure  are 
observed.  However,  the  initial  material  state  is 
not  restored.  That  is,  the  material  has  undergone 
permanent  changes  in  its  molecular  level  topology, 
bond  structure,  and  density.  This  can  be  seen  by 
comparing  results  displayed  in  Figs  4(c)  and  (d). 
These  figures  show  that  while  A-  and  />’- labelled 
silicon  atoms  regain  their  fourfold  coordination 
upon  depressurization  to  ambient  pressure,  the 
C-labelled  silicon  atom  retains  its  fivefold 
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Fig.  5  Typical  pressure  versus  degree-of-compression 
results  obtained  in  the  present  molecular-level 
analysis  of  repeated  pressurization /depressuri¬ 
zation  simulations.  In  each  loading  cycle,  pres¬ 
surization  was  carried  out  to  a  higher  peak 
pressure  followed  by  a  complete  depressuriza¬ 
tion  (i.e.  to  the  atmospheric  pressure) 

coordination,  Figs  4(b)  and  (d).  Thus,  the  average 
coordination  number  of  silicon  atoms  in  the 
relaxed-then- depressurized  state,  Fig.  4(d),  is 
higher  than  that  in  the  initial  configuration,  Fig. 
4(c).  Furthermore,  while  the  aforementioned 
twofold  Si-0  ring  was  broken  upon  depressuriza¬ 
tion,  the  threefold  ring  involving  the  C-,  D-,  and 
Fi-labelled  silicon  atoms  survived  this  process.  In 
addition,  a  new  fourfold  ring  involving  the  Si  atoms 
labelled  F,  G,  H,  and  I  was  formed.  The  threefold 
and  fourfold  rings  were  not  present  in  the  initial 
molecular-level  configuration,  Fig.  4(c). 


4.1.2  Pressure  versus  degree-of-compression 
relation 

An  example  of  the  typical  pressure,  P  versus  degree- 
of-compression,  (.i  results  obtained  in  the  present 
molecular-level  analysis  of  high-pressure  irreversible 
densification  of  glass  is  displayed  in  Fig.  5.  The  results 
in  Fig.  5  show  four  loading  pressurization/ depressur¬ 
ization  cycles.  The  first  cycle  does  not  result  in  any 
irreversible  densification  of  glass  since  the  maximum 
pressure  attained  is  not  high  enough.  The  second  and 
third  cycles  yield  irreversible  densification  and  reveal 
that  this  process  is  associated  with  a  nearly  constant 
ca.  4  GPa  pressure  level.  During  the  past  cycle,  irre¬ 
versible  densification  is  completed  so  that  glass 
behaves  as  a  perfectly  elastic  material  when  subjected 
to  any  further  loading. 


To  summarize,  the  examination  of  the  results  dis¬ 
played  in  Fig.  5  revealed  the  following  three  impor¬ 
tant  findings:  (a)  irreversible  densification  begins  at  a 
pressure  level  of  ~4  GPa  and  proceeds  to  full  densifi¬ 
cation  at  a  nominally  constant  pressure;  (b)  irrevers¬ 
ible  densification  is  associated  with  a  density  increase 
of  ca.  5  per  cent;  and  (c)  the  average  rate  of  change  of 
pressure  with  density  (which  scales  with  the  material 
bulk  modulus)  is  not  significantly  different  (and  will 
be  assumed  equal)  in  the  pre-  and  post- densification 
glass  states. 

It  should  be  noted  that  the  aforementioned  find¬ 
ings  (a)  and  (c)  were  found  not  to  be  very  sensitive  to 
the  random  selection  of  the  initial  molecular-level 
configuration  of  glass.  On  the  other  hand,  the  extent 
of  irreversible  densification  was  found  to  vary  in  a  3-7 
per  cent  range  depending  on  the  choice  of  this  con¬ 
figuration.  This  observation  suggests  that  the  volu¬ 
metric  response  of  the  soda-lime  glass  in  its  pre- 
and  post-irreversible  compaction  states  as  well  as 
the  pressure  at  which  irreversible  compaction  initi¬ 
ates/proceeds  are  functions  of  the  global/assembly 
level  thermodynamic  state  of  the  material  while  the 
extent  of  the  irreversible  densification  is  more  sensi¬ 
tive  to  the  local-details  of  the  material  microstruc- 
tural  state. 

4.1.3  Densification- induced  material 
strengthening 

As  mentioned  earlier,  molecular  level  simple- shear 
tests  were  carried  out  in  order  to  assess  the  extent 
of  irreversible  densification-induced  strengthening 
of  glass.  Towards  that  end,  shearing  of  the  compu¬ 
tational  cell  was  carried  out  in  small  increments 
followed  by  energy  minimization  with  respect  to 
the  atomic  positions.  An  example  of  the  typical 
molecular-level  topology  evolution  accompanying 
these  tests  is  displayed  in  Figs  6(a)  to  (d).  To 
help  prevent  the  computational  crystal  from 
settling  into  a  nearby  metastable  higher-energy 
configuration,  a  10  000-step  300  K  NVT  molecular 
dynamics  (equilibration)  run  was  introduced 
between  the  cell  shearing  and  the  energy  minimi¬ 
zation  steps.  This  procedure  yielded  a  plot  of  the 
potential  energy  increase  (relative  to  that  in  the 
initial  optimized  computational  cell)  versus  the 
shear  angle  of  the  computational  cell.  The  shear 
angle  gradient  of  the  potential  energy  increase 
divided  by  the  constant  unit  cell  volume  is  then 
used  to  define  the  material  shear  strength.  This 
procedure  revealed  that  the  shear  strength  of 
glass  in  the  irreversibly  compacted  state  is  higher 
by  a  factor  of  ca.  1.5  relative  to  the  strength  of  glass 
in  the  initial  state. 
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Fig.  6  An  example  of  molecular-level  topology  evolution  accompanying  simple-shear  mechanical 
tests.  The  tests  were  used  to  assess  the  extent  of  irreversible  densification-induced  strength¬ 
ening  of  glass 


4.1.4  Kinetics/dynamics  of 
irreversible- densification  processes 

As  correctly  pointed  out  by  one  of  the  reviewers  of  this 
manuscript,  the  observed  microstructural  changes, 
which  result  in  irreversible  densification  take  place 
at  room  temperature.  Conventionally,  room  temper¬ 
ature  is  considered  to  be  associated  with  a  low  value 
of  the  material  homologous  temperature  and,  in  the 
case  of  soda-lime  glass,  lies  several  hundreds  of 
degrees  below  the  glass-transition  temperature. 

Hence,  it  is  quite  unusual  that  significant  micro- 
structural  changes  can  occur  in  soda-lime  glass  at 
room  temperature.  However,  it  should  be  recognized 
that,  in  the  present  case,  the  glass  is  subjected  to  high 
levels  of  pressure,  which  greatly  increases  the  poten¬ 
tial  energy  of  this  material  in  its  original  microstruc¬ 
tural  state.  In  other  words,  as  the  pressure  is  increased 
other  microstructural  states  of  glass  become  thermo¬ 
dynamically  preferred.  In  the  lower  range  of 


pressures,  while  the  alternative  microstructural 
states  may  be  thermodynamically  preferred,  the  pres¬ 
ence  of  an  activation-energy  barrier  prevents  (except 
at  sufficiently  high  temperatures)  the  phase  transi¬ 
tion  of  the  material  from  its  initial  microstructural 
state  to  the  thermodynamically  preferred  microstruc¬ 
tural  state.  In  the  upper  range  of  these  pressures,  on 
the  other  hand,  the  activation  energy  barrier  disap¬ 
pears  (at  least,  locally)  so  that  the  phase  transition 
and  the  associated  irreversible  densification  can 
take  place  at  very  high  rates  and  without  any  assis¬ 
tance  from  thermal  activation. 

4.2  JH2  EOS  and  strength-model  modifications 

The  molecular-level  modelling  results  attained  in  the 
previous  section  were  used  to  modify  the  original  JH2 
EOS  and  the  strength  model  (as  presented  in  section 
3).  A  brief  description  of  these  modifications  is  given 
in  the  remainder  of  this  section. 
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4.2.1  Modifications  to  the  EOS 

An  examination  of  the  polynomial  EOS  for  glass  as 
implemented  in  the  ANSYS/Autodyn  material 
library  revealed  that  the  extent  of  non-linearity  is 
quite  small,  and  hence,  the  P  versus  /x  relationship 
was  simplified  using  a  linear  EOS.  In  accordance 
with  the  results  obtained  in  the  previous  section, 
the  following  additional  simplifications/assump¬ 
tions  were  made: 

(a)  the  bulk  modulus  was  assumed  to  take  the  same 
(constant)  value  in  both  the  initial  and  com¬ 
pacted  states  of  glass,  as  well  as  in  all  intermedi¬ 
ate  glass  states; 

(b)  irreversible  densification  of  glass  is  assumed  to 
take  place  at  a  constant  level  of  pressure. 

In  accordance  with  the  molecular-level  computa¬ 
tional  results  presented  in  the  previous  section,  this 
level  of  pressure  was  set  to  a  value  of  4  GPa. 

A  schematic  of  the  changes  made  in  the  P  versus  /x 
relationship  is  given  in  Fig.  7.  The  aforementioned 
modifications  in  the  JH2  EOS  were  implemented  in 
the  ‘ MDEOS_USER_l.f90’  material  user  subroutine 
which  is  then  linked  with  ANSYS/Autodyn  object 
code  to  form  a  new  executable. 

4.2.2  Modifications  of  the  strength  model 

Irreversible  densification  of  glass  is  assumed  to 
increase  its  strength,  both  in  the  damage-free  and 
the  fractured  states.  This  assumption  is  based  on  the 
following  argument:  in  the  case  of  the  damage-free 
glass  state,  irreversible  densification  is  assumed  to 


Fig.  7  A  schematic  the  modifications  made  in  the 
JH2  (over-pressure,  P  versus  degree-of- 
compression,  g)  EOS  in  order  to  account  for 
the  effect  of  high-pressure  irreversible 
densification 


repair  molecular-level  void-like  regions.  In  the  case 
of  the  fractured  glass  state,  irreversible  densification 
is  assumed  to  act  as  a  potent  fragment-deformation 
mechanism  resulting  in  mutually  more-conformal 
fragments.  The  resulting  increase  in  inter-fragment 
contact  surface  areas  would  then  give  rise  to  an 
increase  in  the  friction-controlled  strength  of  fractured 
glass.  An  examination  of  equations  (26)  and  (27) 
revealed  that  the  parameters  A  and  B  control  glass 
strength  in  the  two  states,  respectively.  Following  the 
aforementioned  molecular-level  findings  regarding 
densification-induced  strengthening  of  the  glass, 
parameters  A  and  B  are  increased  by  a  factor  of  1.5 
for  the  fully  densified  glass  material.  At  the  intermedi¬ 
ate  levels  of  glass  densification,  parameters  A  and  B  are 
assumed  to  be  linearly  related  to  the  degree  of  irrevers¬ 
ible  densification.  The  aforementioned  modifications 
in  the  JH2  strength  model  were  implemented  in  the 
‘ MDSTR_USER_l.f90’  material  user  subroutine  which 
is  then  linked  with  ANSYS/Autodyn  object  code  to 
form  a  new  executable. 

4.3  Effect  of  high-pressure  irreversible 
densification  on  ballistic-impact 
resistance  of  glass 

In  the  previous  section,  the  molecular-level  compu¬ 
tational  results  regarding  high-pressure  irreversible 
densification  of  glass  were  used  to  modify  the  EOS 
and  the  strength  relations  within  the  JH2  continuum 


Fig.  8  (a)  A  schematic  of  the  glass  plate  impacted  by  a 

solid  right- circular  cylindrical  projectile  and  (b) 
an  example  of  the  computational  mesh,  based 
on  the  first-order  four-node  finite  elements, 
used  to  model  the  projectile  and  the  glass  plate 
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model  for  this  material.  In  this  section,  a  simple  trans¬ 
verse  impact  of  a  monolithic  glass  plate  by  a  solid 
right- circular  cylinder  (FSP)  is  analysed  computa¬ 
tionally  in  order  to  assess  if  the  aforementioned  mod¬ 
ifications  in  the  JH2  model  have  any  significant  effect 
on  the  temporal  evolution  and  spatial  distribution  of 
damage  within  the  glass  plate  during  impact  and  on 
the  overall  glass-plate  penetration  resistance. 

4.3.1  Problem  formulation 

A  schematic  of  a  circular  disc-shaped  glass-plate 
(25.0  mm  thick  and  75.0  mm  radius)  impacted  by  an 
FSP  (12.7  mm  diameter  and  12.7  mm  length)  problem 
is  provided  in  Fig.  8(a).  Owing  to  the  attendant  sym¬ 
metry  of  the  problem,  an  axisymmetric  2D  formula¬ 
tion  was  employed.  An  example  of  the  finite- element 
mesh  based  on  the  first-order  four-node  finite  ele¬ 
ments  is  displayed  in  Fig.  8(b).  It  is  seen  that  a  finer 
mesh  is  used  in  the  FSP  and  in  the  region  of  the  glass 
plate  which  is  most  severely  affected  during  impact/ 
penetration  of  the  plate  by  the  FSP. 

The  FSP  is  assumed  to  be  made  of  AISI/SAE  4340 
steel  and  this  material  is  modelled  using  a  linear  EOS, 
the  Johnson-Cook  strength  model,  the  Johnson- 
Cook  failure  model,  and  an  erosion  model  based  on 
the  equivalent  geometrical  (i.e.  elastic  +  plastic  + 
crack)  instantaneous  strain.  The  details  regarding  for¬ 
mulation  and  parameterization  of  this  model  can  be 
found  in  our  previous  work  [28] . 

All  the  calculations  carried  out  in  this  study  were 
done  using  AN SYS/ Autodyn,  a  general  purpose  non¬ 
linear  dynamics  modelling  and  simulation  software 
[26] .  A  transient  non-linear  dynamics  problem  is  ana¬ 
lysed  within  AN SYS /Autodyn  by  solving  simulta¬ 
neously  the  governing  partial  differential  equations 
for  the  conservation  of  momentum,  mass  and 
energy  along  with  the  material-model  equations 
and  the  equations  defining  the  initial  and  the  bound¬ 
ary  conditions.  The  equations  mentioned  above  are 
solved  numerically  using  a  second- order  accurate 
explicit  scheme  and  one  of  the  two  (Lagrange  or 
Euler)  basic  mathematical  approaches/processors. 
The  key  difference  between  the  two  processors  is 
that  within  the  Lagrange  processor  the  numerical 
grid  is  attached  to  and  moves  along  with  the  material 
during  calculation  while  within  the  Euler  processor, 
the  numerical  grid  is  fixed  in  space  and  the  material 
moves  through  it.  In  this  study,  both  the  FSP  and  the 
glass  plate  are  modelled  using  the  Lagrange  proces¬ 
sor.  In  our  recent  work  [28],  a  brief  discussion  was 
given  of  how  the  governing  differential  equations 
and  the  material-model  relations  define  a  self- 
consistent  system  of  equations  for  the  dependent 
variables  (nodal  displacements,  nodal  velocities, 


material- element  densities,  and  material-element 
internal  energy  densities).  In  the  same  reference,  a 
brief  overview  was  provided  of  the  contact/interac- 
tion  and  erosion  algorithms  implemented  in  AN SYS/ 
Autodyn. 

The  initial  conditions  are  defined  is  such  a  way  that 
the  FSP  is  assigned  a  downward  velocity  while  the 
glass  plate  is  assumed  to  be  stationary.  To  reduce  the 
effect  of  reflection  of  the  shock  waves  at  the  (outer) 
hoop  surface  of  the  glass  plate,  the  so-called  ‘transmit’ 
boundary  conditions  were  applied  to  all  the  nodes 
residing  on  this  surface.  The  transmit  boundary  con¬ 
ditions  enable  propagation  of  the  stress  waves  across 
the  hoop  surface  without  reflection  [28] . 

4.3.2  Glass-plate  impact  computational  results 

An  example  of  the  results  pertaining  to  the  temporal 
evolution  and  spatial  distribution  of  deformation  and 
damage  within  the  glass-plate  and  the  FSP  is  dis¬ 
played  in  Figs  9(a)  to  (d).  Various  deformation/failure 
models  are  colour  coded  to  improve  clarity  of  the 
results.  The  results  displayed  in  Figs  9(a)  to  (d)  show 
that  shear-induced  damage /failure  is  found  mainly  in 
the  region  underneath  the  FSP,  tensile  failure  occurs 
at  the  glass-plate  back-face  (and  is  caused  by  reflec¬ 
tion  of  the  compressive  stress  wave  at  this  face),  and 
the  FSP  undergoes  extensive  plastic  deformation  (but 
little  or  no  failure) .  These  findings  are  fairly  common 
and,  hence,  will  not  be  discussed  any  further.  It 
should  be  noted,  however,  that  the  results  displayed 
in  Figs  9(a)  to  (d)  pertain  to  an  FSP  initial  velocity  of 
600  m/s  and  the  original  (unmodified)  JH2  glass- 
material  model. 

When  the  JH2  material  model  was  modified  to 
include  the  effect  of  high-pressure  irreversible  densi- 
fication  and  the  FSP  initial  velocity  was  maintained  at 
600  m/s,  an  identical  temporal  evolution  and  spatial 
distribution  of  the  material  deformation  and  failure 
was  found  as  those  in  Figs  9  (a)  to  (d) .  This  finding  was 
readily  rationalized  by  establishing  that  the  pres¬ 
sure  within  the  glass-plate  never  reached  the  level 
of  4  GPa  required  for  the  irreversible  densification  to 
commence. 

When  the  FSP  velocity  was  increased  to  1000  m/s, 
the  effect  of  irreversible  densification  was  observed. 
An  example  of  this  effect  is  given  in  Figs  10(a)  and  (b), 
in  which  distribution  of  deformation  and  damage 
within  the  glass-plate  and  the  FSP  are  given  for  the 
cases  of  original  and  modified  JH2  glass  material 
models,  respectively.  The  results  displayed  in  these 
figures  show  that  the  extent  of  damage  in  the  FSP 
increased,  while  that  of  damage  in  the  glass-plate  is 
reduced  for  the  case  of  the  modified  JH2  material 
model.  Also,  the  FSP  exit  velocity  was  found  to  be 
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Fig.  9  An  example  of  the  results  pertaining  to  the  temporal  evolution  and  spatial  distribution  of 
deformation  and  damage  within  the  glass-plate  and  the  FSP 


reduced  from  469.1  to  460.2  m/s  as  a  result  of  glass 
material  model  modifications.  This  is  a  relatively 
small  decrease  (ca.  2  per  cent)  in  the  FSP  residual 
velocity  which  may  suggest  that  the  irreversible  com¬ 
paction  of  glass  may  not  be  a  potent  ballistic-resis¬ 
tance  enhancing  mechanism  in  glass. 

To  further  explore  the  effect  of  irreversible  densifi- 
cation,  a  parametric  study  was  set  up  in  which  the 
key  parameters  of  this  process  were  varied  within 
physically  realistic  limits.  Specifically,  the  pressure 
at  which  densification  starts  was  varied  in  a  range 
1-10  GPa  range,  density  increase  between  5  and 
15  per  cent  and  the  densification-induced  strength¬ 
ening  extent  (i.e.  the  multiplier  of  A  and  B  parame¬ 
ters)  in  a  1. 0-2.0  range.  This  study  revealed  that, 
at  an  FSP  velocity  of  1000  m/s,  the  FSP  exit  velocity 
can  be  reduced  from  469.1  to  346.4  m/s  (a  26  per  cent 
reduction).  The  optimal  set  of  irreversible- 
densification  parameters  was  identified  as:  densifica¬ 
tion  pressure  1  GPa,  density  increase  5  per  cent,  and  a 


strengthening  extent  of  2.0.  Figure  10(c)  reveals  spa¬ 
tial  distribution  of  deformation/ damage  for  the  mod¬ 
ified  JH2  glass  model  and  the  optimal  set  of 
irreversible -densification  parameters.  It  is  seen  that 
there  is  a  significant  reduction  in  damage  to  the  glass 
material  directly  under  the  projectile  due  to  the 
increased  strength  of  the  densified  glass.  This  reduc¬ 
tion  in  failed  material  under  the  projectile  promotes 
an  increase  in  the  extent  of  plastic  deformation  expe¬ 
rienced  by  the  projectile  which  leads  to  spreading  of 
its  material  over  a  larger  area.  These  observed  effects 
combine  to  allow  for  a  significant  decrease  in  the 
kinetic  energy  of  the  projectile  with  respect  to  the 
non-optimized  glass  material  case. 

Additional  simulations  were  next  carried  out  in 
order  to  establish  if  irreversible-compaction  of  glass 
(if  associated  with  the  aforementioned  optimal  set  of 
parameters)  can  improve  the  ballistic-resistance  of 
glass  at  lower  FSP  velocities.  It  was  found  that  at  an 
FSP  initial  velocity  of  600  m/s,  the  FSP  residual 
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Fig.  10  The  effect  of  the  JH2  material-model  modifi¬ 
cations  on  the  spatial  distribution  of  damage 
within  a  soda-lime  glass  plate  impacted  with  a 
12.7-mm  diameter/ 12. 7-mm  height  solid 
right-circular  cylinder  at  an  initial  velocity  of 
1000  m/s:  (a)  the  original  model;  (b)  the  model 
modified  using  molecular-level  finding;  and 
(c)  the  modified  using  an  optimal  set  of  high- 
pressure  irreversible-densification  parameters 


velocity  was  reduced  from  240.4  to  154.2  m / s.  Thus,  if 
the  irreversible-densification  parameters  of  glass  can 
be  set  to  their  optimal  values,  then  the  ballistic-pene¬ 
tration  resistance  can  be  improved  for  a  range  of  FSP 
impact  velocities. 

4.3.3  Experimental  validation 

While  there  are  numerous  reports  of  experimental 
investigations  aimed  at  assessing  the  ballistic  perfor¬ 
mance  of  soda-lime  glass  [12-15,  31],  there  are  no 
open-literature  reports  regarding  the  operation  or 
the  role  of  irreversible  densification  in  improving 
the  ballistic  performance  of  soda-lime  glass.  To 
address  this  problem,  a  companion  experimental 
investigation  is  underway  at  the  Army  Research 
Laboratory,  Aberdeen,  Proving  Ground,  Maryland. 
Preliminary  findings  obtained  in  the  post-mortem 
X-ray  diffraction  studies  of  glass  fragments,  carried 
out  as  part  of  this  investigation,  show  the  presence 
of  a  higher  density  soda-lime  glass  microstructural 
state  and  can  be  found  in  reference  [32] . 


Table  1  Johnson-Holmquist  2  material  model  param¬ 
eters  for  soda-lime  glass 


Parameter 

Symbol 

Unit 

Value 

EOS:  linear 

Bulk  modulus 

K, 

Pa 

4.54el0 

Strength  model:  Johnson-Holmquist  2 

Shear  modulus  G 

Pa 

3.04el0 

HEL 

&HEL 

Pa 

5.95e9 

Intact  strength  constant 

A 

N/A 

0.93 

Intact  strength  exponent 

N 

N/A 

0.77 

Strain  rate  constant 

C 

N/A 

0.003 

Fractured  strength  constant 

B 

N/A 

0.35 

Fractured  strength  exponent 

M 

N/A 

0.40 

Max  fracture  strength  ratio 

- 

N/A 

0.50 

Failure  model:  Johnson-Holmquist  2 

Hydro  tensile  limit  T 

Pa 

-3.5e7 

Damage  constant  1 

A 

N/A 

0.053 

Damage  constant  2 

A 

N/A 

0.85 

Bulking  constant 

P 

N/A 

1.0 

4.3.4  Materials -by -design  approach 

In  summary,  the  results  obtained  in  this  study  show 
that  high-pressure  irreversible  densification  of  glass 
can  become  a  potent  ballistic-resistance  enhancing 
mechanism  provided  (via  chemical  modifications, 
and  various  thermo-mechanical  and  chemical  treat¬ 
ments)  the  densification  parameters  can  be  set  to  or 
near  optimal  values.  In  our  ongoing  work,  molecular- 
level  modelling  is  being  extensively  used  to  help 
better  define  such  glass-modifying  strategies/proce¬ 
dures.  Specifically,  within  the  molecular-level  model¬ 
ling  framework,  the  effect  of  the  glass  chemistry 
and  its  random-network  structure  (as  characterized 
by  a  set  of  random-network  microstructural  parame¬ 
ters  [29,  30])  on  the  irreversible-densification  charac¬ 
teristics  (i.e.  the  onset  pressure  and  the  extent  of 
density  change)  is  being  examined.  In  addition,  a 
link  is  being  established  between  the  molecular- 
level  modelling  parameters  and  the  JH2  continuum- 
level  material  parameters.  Table  1  summarizes  the 
JH2  model  parameters  and  lists  their  typical  values 
for  the  case  of  a  ‘non- transforming’  soda-lime  glass. 
In  this  study,  EOS  is  modified  to  account  for  the 
irreversible  densification  transition  and  two  strength 
parameters  (parameters  A  and  B)  were  modified. 
In  our  on-going  work,  an  attempt  is  being  made  to 
identify  additional  JH2  model  parameters  which 
require  modifications  as  the  soda-lime  glass  chemis¬ 
try  and  the  random-network  structure  are  being  mod¬ 
ified.  It  is  hoped  that  this  will  become  a  fruitful 
exercise  of  the  so-called  ‘materials-by-design’  con¬ 
cept  in  which  component-level  experimental  and 
computational  investigations  are  carried  out  in 
order  to  identify  the  envelope  of  optimal  material 
properties.  This  is  then  followed  by  an  extensive 
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experimental/computational  procedure  aimed  at 
designing  and  synthesizing  the  materials  with  the  tar¬ 
geted  set  of  properties/performance  attributes. 

5  CONCLUSIONS 

Based  on  the  results  obtained  in  this  study,  the  fol¬ 
lowing  main  conclusions  can  be  made. 

1.  Molecular-level  modelling  of  soda-lime  glass 
revealed  the  occurrence  of  an  irreversible-densifi- 
cation  process  when  the  pressure  exceeds  ca. 
4  GPa.  Close  examination  of  molecular-level  topol¬ 
ogy  revealed  that  this  process  is  associated  with  an 
increase  in  the  average  coordination  number  of  the 
silicon  atoms,  and  the  creation  of  two  to  fourfold 
(smaller,  high  packing-density)  Si-O  rings. 

2.  Modifications  of  the  continuum-level  material 
model  for  glass  to  include  the  effect  of  irreversible 
densification  resulted  in  minor  improvements  in 
the  ballistic-penetration  resistance  of  glass  and 
only  for  high  projectile  initial  velocities.  Also,  it  is 
observed  that  the  main  threat  to  a  soda-lime  glass- 
based  transparent-armour  panel  during  the  ballis¬ 
tic  impact  arises  from  the  normal/tensile  and  less 
from  deviatoric/ shear  stresses. 

3.  A  parametric  study  involving  variations  of  the  key 
irreversible -densification  parameters  within  physi¬ 
cally  realistic  limits  revealed  that  the  optimal  com¬ 
bination  of  these  parameters  can  result  in 
substantial  improvements  in  the  ballistic-resistance 
of  glass  over  a  wide  range  of  projectile  velocities. 

4.  It  is  suggested  that  various  chemical-modification 
and  thermo-mechanical  treatment  strategies 
should  be  employed  in  order  to  attain  this  optimal 
set  of  irreversible- densification  parameters.  This  is 
an  example  of  the  materials-by-design  concept 
within  which  an  optimal  combination  of  material 
parameters  is  identified  in  order  to  maximize  com¬ 
ponent-level  performance. 
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